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We show that two identical solid-state qubits can be made 
fully entangled (starting from completely mixed state) with 
probability 1/4 just measuring them by a detector, equally 
coupled to the qubits. This happens in the case of repeated 
strong (projective) measurements as well as in a more realistic 
case of weak continuous measurement. In the latter case the 
entangled state can be identified by a flat spectrum of the 
detector shot noise, while the non-entangled state (probability 
3/4) leads to a spectral peak at the Rabi frequency with the 
maximum peak-to-pedestal ratio of 32/3. 

Prospective solid-state realizations of quantum 
computerail have many potential advantages over other 
proposed realizations, including more natural scalabil- 
ity, simple electrical control of parameters, and use of 
very well developed technology. A number of theoretical 
proposals on solid-state quaxitum computers have been 
put forward (see, e.g. RefsBlj) and spectacular experi- 
mental results have begn achieved, including demonstra- 
tions of chaxge qubitsQ using single-Cooper-pair boxes, 
flux qubitsEfcl using superconducting loops interrupted by, 
Josephson junctions, and combined charge-flux qubitst2l 
(the deaionstrated quality factor was as high as 25,000 
in RefJ13). Obviously, next important experimental step 
is the demonstration of entangled solid-state qubits. 

The qubits can be made entangled by their direct in- 
teraction if the interaction depends on the states of the 
qubits and at least one qubit is in a coherent superpo- 
sition before the interaction. In this letter we will dis- 
cuss an alternative way, when two solid-state qubits are 
made entangled just by their simultaneous measurement 
by one detector, which thus provides an indirect coupling 
between qubits. The procedure works with a probability 
less than unity, and in this respect it is some'wdiat similar 
to the operation of conditional quantum gatestil based on 
simple linear optical elements. Our procedipe can also 
be thought of as an entanglement purificatioiJla, in which 
instead of trading the number of qubit pairs for better 
entanglement, we trade the success probability while no 
initial entanglement is needed. One more analogy is the 
preparation of entangled atpsns in an optical cavity by 
monitoring the cavity decayli^. 

In contrast to qubits represented by photons, which are 
physically destroyed by the acts of measurement, solid- 
state qubits only change their state due to measurement, 
that allows somewhat more freedom in designing quan- 
tum operations. On the other hand, it is quite difhcult 
to realize simple projective measurements of solid-state 
qubits because of typically weak coupling with detector. 
Therefore, instead of simple abrupt collapse, we have 



to deal with dephasing^like processes in the case of en- 
semble measuj^naeptaij or with the continuous (weak) 
measurementaljiia in the case of single qubits. 

The theory of non-averagftd, continuous measurement 
of single solid-state qubitalJ'Ea is under active develop- 
ment during last 4 years and has already produced a 
number|-e|f experimental predictions-/ for some of them, 
see Refli^). One of the predictionsc3 is related to a di- 
rect continuous measurement of a qubit Rabi oscillations 
by a weakly coupled detector. The theoretical result is 
that the qubit oscillations can be evidenced by the peak 
of the detector current spectral density at the Rabi fre- 
quency; however, the peak height cannot be larger than 
4 times the noise pedestal pt|his fact seems to have recent 
experimental confirmationEj). 

In this letter we consider the case when two identical 
qubits perform Rabi oscillations, which are continuously 
measured by a detector equally coupled to two qubits. 
We will show that the system is gradually collapsed into 
one of the two regimes: either qubits become fully entan- 
gled (Bell state) that can be identified by absence of the 
spectral peak of the detector current, or qubits state fall 
into the orthogonal subspace that can be identified by the 
Rabi spectral peak, which for an ideal detector is 32/3 
times higher than the noise pedestal. The probabilities 
of two scenarios are 1/4 and 3/4, respectively (starting 
from completely mixed state) , so on average the peak-to- 
pedestal ratio is equal to 8, twice as large as for a single 
qubit. 

Figure |l| shows possible realizations of our setup. In 
the first realization [Fia. E|(a)] each qubit is made of 
a double quantum dotc3 (DQD), occupied by a single 
electron, while the detector is a quantum point con- 
tact (QPC) located in between DQDs. The second pos- 
sible realization [Fig. 0(b)] is, based on single-Cooper- 
pair-boxes (SCPB) as qubitsQ, which are measured by 
a single-electron transistor (SET). Other possibleureal- 
izations (not shown) can be_hased on flux qubitstltl or 
combined charge-flux qubitalJ. 

In the Hamiltonian of the system, 

Ti- = Ti.QB +'Hdet + 'HiNT, (1) 

the first term describes two qubits alone. Tigs = 
(ea/2)(a|ax - ajat) -t- Ha{a^^ai + ajat) 4- (eh/2)(6j6x - 

b]-b^) + Hi,{b^^bi + b^^b{), where Sa and Sb are energy asym- 
metries, which are neglected in this letter, Sa = Sb = 0, 
the amplitudes Ha and Ht describe the tunneling within 
qubits (we assume Ha = Hf, in the most of this letter), 
and the direct interaction term J7a|a|6|6| is neglected. 
The frequencies of free Rabi oscillations of qubits, Qa — 



1 




V V, 



'gb (b) 



DQDfl QPC DQDfe 



SQPBa SET SQPBfo 



FIG. 1. Schematic of two qubits measured by an equally 
coupled detector, (a): Realization based on double quan- 
tum dots measured by a quantum point contact, (b): real- 
ization based on single-Cooper-pair-boxes measured by a sin- 
gle-electron transistor. Measurement can entangle qubits. 
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(we use ?i = 1) obviously coincide, f2a = = 17 
[f2 = (Oq -I- rJfc)/2] if Hcr-T^ Hb- For simplicity wc 
limit ourselves by the casecj of DQD qubits, measured 
by a low transparency QPC (though generalization to 
other cases is simple), so that the detector Hamiltonian 

is HdET = J2l Eicjci + J2r Erclcr + J^l.r Ticjcr + cjc/) 

and the interaction term is "Hint = X); r ^-^^('^t^T ^ 
aja|)(c|cr + 4c/) + E;,r Arf,(5|5| - 6|5x)(c|'cr + 4q); 
we will be mostly interested in the case of equal coupling, 
AT, = ATfc. 

The four basis states of two qubits, |1) ^jtatb), 
|2) =|Taifc), |3) =UaT6), |4) =|iai6), correspond to 
4 values of the average current through the detector: 
h.2,3A = 2tt{T ± ATa ± ATbfpiprC^V, where V is the 
QPC voltage and pi(j.) are densities of states. It is impor- 
tant that in the case of equal coupling two currents co- 
incide, I2 = I3 = I23, so the measurement cannot distin- 
guish between states |2) and |3). Besides the "measure- 
ment" basis, it is convenient to introduce also the Bell ba- 

sis: = (ITaib)- \lah))/V2, |2)^ = i\UU)- \laib 

))/^/2, |3)^ ^ (|Tai6>+ \iaU))/V2, and |4)^ ^ {\UU 
)+ |iai&))/\/2. Notice that and |2)^ are eigenstates 
of Hqb if Ha = Hb, while the states |3)^ and |4)-^ are 
transformed by Hqb as cos(fit-F0) 1 3)^ -i sin(f7t-h0) |4) ^ . 

Before considering continuous measurements, let us 
discuss a simpler case of a sequence of "orthodox" projec- 
tive measurements which can be realized if the coupling 
with the detector is strong {C :s> 1, see below) and detec- 
tor voltage is applied during short time intervals. Since 
the states |2) and |3) are mutually indistinguishable, the 
the two-qubit density matrix p is projected each time into 
one of three subspaces, corresponding to states |23), 
and 1 4) (we use notation |23) for the subspace spanned 
by 1 2) and |3)). The sequence of measurements separated 
by time periods At can be described by such projections 
separated by intervals of unitary evolution due to Hqe- 

Suppose the first measurement resulted in the cur- 
rent I23, then the system is projected into |23) subspace, 
which is also a subspace |13)^ in the Bell basis. If the 
state would be exactly |1)^ (which does not evolve un- 
der Hqb), then all subsequent measurements would give 



the same result I23 and the state |1)^ would remain un- 
changed. However, if the two qubits would be in the state 
1 3)^, then the next measurement would result again in 
I23 only with probability p = (cos fiAt)^, while the prob- 
abilities of results Ii and I4 would be (1 — p)/2 each. 
Therefore, if a sufficiently long sequence of current mea- 
surements repeatedly gives the result I23, the two-qubit 
density matrix p purifies and becomes close to the fully 
entangled state 

Simple analysis shows that after N successful measure- 
ments (all results are 123) 



pfi(O) 



Pfl(0)+pf3(0)(cOsl7At)2(A'-l) 



(2) 



where pfi(O), p^(0), and pfi{N) are the correspond- 
ing density matrix elements in the Bell basis before and 
after the sequence of measurements, while the prob- 
abihty of a successful sequence is P{N) — pfi(O) -t- 
/9f3(0)(cosriAt)2(^-i). For large N the difference be- 
tween the obtained state and state jl)"^ becomes expo- 
nentially small, while the probability of success is close 
to pfi{0), which is equal to 1/4 for the fully mixed initial 
state pij^mix = pfj^mix — ^ij I ^ (this state is a direct prod- 
uct of fully mixed states of each qubit) . The purification 
rate depends on At, and is the fastest when At is close 
to (1/4 -I- fc/2)27r/r2 (fc is an integer), which is a regipafi 
opposite to the quantum nondemolition measurementscj. 

Notice that if some measurement in a sequence resulted 
in the current 1\ or I 4, then the subsequent measure- 
ments (for general At) can result in any current l\, I^, 
or I23; however, pf^ remains exactly zero (therefore, long 
sequences of I23 results become extremely improbable). 
Hence, to obtain the Bell state |1)'^, one have to apply 
some perturbation which mixes two subspaces (for exam- 
ple, apply a noise which affects Eq and/or Sb) and repeat 
the procedure. Assuming a well-mixed case, the probabil- 
ity 1 — (3/4)^^ to obtain the state |1)^ can be made arbi- 
trary close to unity by allowing sufficiently large number 
M of attempts. 

The procedure can obviously be used for the prepara- 
tion of entangled states in a solid-state quantum com- 
puter, so it is important to discuss what happens if the 
conditions Ha = Hb and I2 = I3 are not satisfied ex- 
actly. In the case of slightly different Ha and Hb, Eq. 
(^ changes insignificantly (cosfiAt should be replaced 
with cosr2At/cos2~iArJAt, where Afi ee fj, - rjf,), 
however, the probability of an A'^-long successful se- 
quence becomes P{N) = pfi(0)(cos 2-iAf7At)2(^-i) + 
p^(0)(cos ilAt)^*^^"^^ and decreases to zero at ^ cx3. 
Estimating the average length of a successful sequence, 
N ~ (sin 2~^AJ7At)~^, one can estimate a typical devi- 
ation from state 11)^, 



(cosl7At/cos2-iA17At)2/(«'"2-iAnAt)^^ 

which is as small as ^ exp[— (2f7/Ar2)^] if At ^ 
(Quantum Zeno regime) and even smaller, ~ 
(cosl7At)(32A')(^VAi2)=^^ if is close to Tr/2n. 
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To analyze the effect of a small difference between I2 
and I3 because of slightly different coupling, we have 
to use the standard theoryll3 of weak quantum mea- 
surements and take into account the detector shot noise 
Si = 2eli. We assume that during short measurement 
interval St the currents Ii and £4 can be unambiguously 
identified, while current I2 and I3 are almost indistin- 
guishable and small corresponding signal-to-noise ratio 
is characterized by the parameter e = {I2 — h)^ /4:D <C 1 
where D = 823/261 is the variance of the measured 
noisy current. Each successful measurement tends to 
shift the state towards either |2) or |3) and so decreases 
the amount of entangled state that competes with 

the purification due to Eq. (||) and leads to iterative for- 
mula pfj^{N+l) ~ pfi(A^)-e/4-K[l-pfi(iV)](sinr2Ai)2 
valid when pf^ is close to unity (exact formula is longer). 
Therefore, a typical deviation from pure entanglement 
1 — pfi ~ e/4(sinriAt)^ scales as e. 

Now instead of instantaneous measurements let us 
consider a more realistic case of a continuous measure- 
ment, realized when the detector voltage is applied all 
the time.|— |For the analysis we will use the Bayesian 
formalismliJ assuming weakly responding linear detect- 
ing regime, |A/a,6| < Ii, Ala = h - h = h - h, A/& = 
Ii~ I2 — I3 ~ Ii, and concentrating on the case of sym- 
metric weak couphng, Ca « Cb < 1, Ca = {AlaY / S^Ha, 
Cb = (A/f,)^/S'oi?6, where the frequency-independent de- 
tector shot noise spectral density Sq does not depend 
significantly on the qubits state. 

The evolution of the two-aubit density matrix p can 
be described by the equations (in Ito representation) 

(lit) - J2k Pkklk){lt + - 2 P'^^^^'' ^ 
-((/, - /,)V45o + 7*,) P^o ~ if^QB, (4) 

where the extra dephasing rate 7ij = (?7~^ — ^){Ii — 
IjY/'iSo depends on detector ideality pf^i, (0 < ?y < 1) 
and vanishes for the QPC as a detectortZl (77 — 1); how- 
ever, this term is important, for example, for the SET. 
To simulate individual realizations of the random mea- 
surement process, the noisy detector current I[t) can be 
calculated as 

i{t) = m +Y.kp^^^^^ (5) 

where £^{t) is a white noise with spectral density 5*0. No- 
tice that the averaged dynamics (master equation) can 
be obtained by averaging over noise (^{t) that eliminates 
the first term in Eq. (0). 

We have, performed extensive Monte Carlo 
simulationgl3 and found the following (Fig. H). In the 
symmetric case. Ha = Hb, Ca = Cb — C [C = {Ca +Cb)/2, 
we mostly used C between 1 /4 and 1] , any initial state ei- 
ther evolves eventually into the fully entangled Bell state 
|1)^ (Pn ~^ 1) or ends up in the orthogonal subspace 
{pf^ — > 0) performing oscillations within this subspace 
so that the "signal" z = pu — P44 — 2Rep^ (which af- 
fects the detector current) oscillates with frequency ft 
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FIG. 2. Two numerical realizations of pfi evolution start- 
ing from the fully mixed state. The upper (solid) line illus- 
trates the scenario of collapse into fully entangled Bell state 
1)^, while the lower (dotted) line shows a collapse into the or- 
thogonal subspace. Two insets show the corresponding spec- 
tral densities Si{uj) of the detector noise (solid/dashed lines 
are the numerical/analytical resuls). 

and fluctuating amplitude within the range from to 1. 
Within the accuracy consistent with the number of tri- 
als, the probability of evolving into state coincides 
with pf^{0), similar to the case of orthodox measurement 
sequence. For 77 = 1 the two-qubit state eventually be- 
comes pure independentlv of an initial choice (similar 
to the one-qubit casellZl'ES) . This fact is obvious for the 
state while in the oscillating scenario the surviv- 

ing nondiagonal matrix elements in the Bell basis satisfy 
equations (Repi^f = P22P44, (Impfa)^ = pi2P33^ and 
(Imp;^)^ = pfjPii- In the state |1)^ the numerically 
calculated spectral density of the detector current is flat 
and equal to Sq (the signal z is zero), while in the oscil- 
lating state it exhibits a peak (lower inset in Fig. |2|) at 
frequency with the peak height close to the analytical 
result (32/3)S'o (see below). 

The fact of collapsing eventually either into the state 
or into the orthogonal subspace can be understood 
using an analogy with the sequential measurement case 
considered above, and is caused by the fact that neither 
unitary evolution due to Tigs nor nonunitary evolution 
due to measurement mixes two subspaces [see Eq. (^]. 
The probability of two scenarios should be equal to the 
contribution of two subspaces to the initial density ma- 
trix, pfi(O) and l-pfi(O), since ensemble averaged value 
(pfi(O) does not change with time (as follows from the 
master equation). 

To find analytically the spectral density of the detec- 
tor current-^or the oscillating state, we have used two 
approacheg23 leading to the same result. The first one is 
based on the master equation and collapse ansatz. Us- 
ing classical equation /(<) = zAI + £,{t), we calculate 
the current correlation function Ki(t) = (/(0)/(t)) as 
Ki{t > 0) = {AI)^Ki{T), while Ki{T) is calculated in 
the following way. At time r = the two-qubit state is 
collapsed into one of the three basis states of the sub- 
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space: |1) (corresponding to Ji and z = 1), |4) (corre- 
sponding to /4 and z = —1), or |3)'^ (corresponding to 
J23 and 2: = 0). The probabilities of these coUapses are 
1/3 each, since for pf^ — the stationary solution of 
the master equation is p\\ = P44 = pf^ = 1/3 (the sta- 
tionary solution is diagonal in the Bell basis as welLas 
in the basis {|1), |4), |3)'^} and has zero entanglements). 
In each of 3 cases, the value of z at time r is obtained 
from the solution of the master equation for two relevant 
components: 



dz/dr — —fly, dy/dr — ilz — Fy, 



(6) 



where y = 2Impf3 and T = 7/-i(A/)V45o. So, z(t) = 
±G'(t), G(r) = exp(-rT/2)[cosiir-|-(r/2fi)siniir] [here 
O = (rj2 - rV4)^/^] in the two first cases while z{t) = 
in the third one. Summing 3 contributions to (z(O)z(r)) 
with probability weights 1/3 each, we obtain — 
(2/3)G(r) and the current spectral density 



5*7(0;) = 5o + 



f]2(A/)2r 
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In the case F ^ 51 the spectral peak at the Rabi fre- 
quency n corresppiids to the Q-factor of Sry/C (similar 
to one-qubit caseE3) and has the peak height equal to 
(32/3)r7S'o, confirming the numerical result (see lower in- 
set in Fig. |). 

The second method of Si{uj) calculation is based on 
the Bayesian Eq. (^ assuming ij = 1 and random evolu- 
tion of a pure state with z = A{t)cos{ilt ^^{t)) [then 
y = A{t)sm{Qt + ^{t))]. In this methoclHd the corre- 
lation between noise ^(0) and evolution of the density 
matrix at later time should be taken into account, so 
Ki{t > 0) = {AI)^K,{t) + AIK^^{t), while correla- 
tion functions Kz{t) and K^z{t) should be calculated 
by averaging of iSulong individual realization over time. 
We have provedEZi that the result for Ki{t) calculated 
by this method coincides with the result of the previous 
method for arbitrary coupling C; however, the formalism 
is much simpler for weak coupling, C ^ 1. In this case 
the stochastic differential equations for A{t) and <I>(t) 
can be averaged over oscillations with frequency apd 
the correlation functions can be calculated analyticallytZl: 
K^{t) = (5/12)G(r) and K^,{t > 0) = G(r)A//4. This 
gives us a natural partition of the relative spectral peak 
height 32/3 into two costributions: "classical" part 20/3 
comes from oscillation^ of the signal z, while the "quan- 
tum" contribution equal to 4 is due to partial collapse 
of p correlated with the detector noise. Comparing this 
partition witit,the partition 4 = 2 + 2 for a one-qubit 
measurementciJ, we see that the classical part grows faster 
than the quantum part when the number of qubits is in- 
creased. 

Numerical simulations show that if the two Rabi fre- 
quencies fla and fib are slightly different, or small dif- 
ference between Ca and Cb is due to asymmetry of cou- 
pling (different Ala and Alb), than the two-qubit den- 
sity matrix p makes rare abrupt jumps between a state 



very close to |1)'^ and the oscillating state. To find 
the switching rate analytically, we have used the master 
equation starting from entangled initial condition pfi = 1 
and calculated the linear term in pf^ (t) dependence at 
t ^ r^i_,(but when pfi(i) is still close to unity). In 



this way§ we have obtained the rate Tb^o = (AfJ) /2F 
of switching from the Bell state to the oscillating state 
due to slightly different Rabi frequencies, and the rate 
Tb^o = {AC/CfT/S when fla^flb, but couplings Ala 
and Alb are slightly different. To find the rate of the re- 
verse switching, notice that the stationary master equa- 
tion has the solution pf^ — Pij^st = ^ij/^ (two subspaces 
are mixed now), therefore the system should spend on 
the average 1/4 of the time in the state |1)^, and so 
Fo— >B = Tb^o /3- The numerical histograms of switch- 
ing time distributions confirm these formulas. Taking 
into account rare switching events, the average spectral 
density of the detector current is given by Eq. (^ mul- 
tiplied by 3/4, so the spectral peak height is equal to 
877^0. 

Finally, we have studied the effect of environmental de- 
phasing, modeling it with two small dcphasing rates 7^ 
and 7b acting separately onto two qubits. This leads 
to a slightly mixed— j9 even for an ideal detector and 
to switching event£3 (similar to the case above) with 
Tb^O = 3Fo^_B = (7a + lb)/'2'- Notice that a control- 
lable weak external noise can be used in a simple feedback 
protocol to restore the entangled state after undesirable 
switching to the oscillating state. 

In conclusion, we have found that the continuous mea- 
surement of two identical solid-state qubits by equally 
coupled detector leads to either full spontaneous entan- 
glement of qubits (Bell state) or to collapse into orthogo- 
nal oscillating state. In the latter case the noisy detector 
current has a spectral peak at the Rabi frequency, while 
in the former case the spectrum is flat. Slight asymmetry 
of the two qubit configuration as well as environmental 
dephasing leads to switching between two regimes. It is 
important to mention that for an experimental observa- 
tion of the phenomenon the quantum ideality rj of the de- 
tector should not necessarily be close to unity; it should 
only be large enough to allow distinguishing the Rabi 
spectral peak with the peak-to-pedestal ratio of 3277/3. 

The work was supported by NSA and ARDA under 
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